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Stefan-Maxwell  model  (SMM)  and  simple  Fields  model  (FM)  type  of  relations  both  including  Knudsen 
diffusion  for  the  calculation  of  species  mole  fraction  distribution  inside  the  porous  anode  of  a  solid  oxide 
fuel  cell  (SOFC)  were  compared  and  it  was  found  that  at  low  current  densities  the  models  agree  well 
but  as  current  increases  the  differences  also  increase.  Based  on  the  findings  an  empirical  correction  is 
proposed  for  the  effective  diffusivity  used  in  Fields  model.  The  corrected  diffusivity  coefficient  gave  better 
agreement  with  the  Stefan-Maxwell  model  and  even  at  higher  current  densities  the  error  is  less  than  5%. 
This  correction  was  implemented  via  a  three-dimensional,  in-house  SOFC  simulation  code  (DREAMSOFC) 
which  uses  Fields  model  type  relations  for  diffusion  flux  calculations.  The  code  also  takes  into  account 
methane  steam  reforming  (MSR)  and  water  gas  shift  ( WGS)  reactions  and  the  electrochemical  oxidation  of 
both  H2  and  CO.  As  an  application,  a  SOFC  button  cell  which  is  being  tested  at  West  Virginia  University  was 
simulated.  The  results  with  and  without  the  proposed  correction  for  effective  diffusivity  are  compared. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  perceived  to  be  the  most  viable 
candidates  for  future  power  plants  due  to  their  high  efficiency  and 
fuel  flexibility.  The  performance  of  a  SOFC  at  high  operating  current 
densities  is  reduced  mainly  due  to  the  mass  transport  losses,  that  is, 
the  demand  for  reactants  exceeds  their  capacity  to  diffuse  through 
the  porous  anode  to  the  reaction  site  at  the  anode-electrolyte  inter¬ 
face.  The  multi-component  mass  transport  in  SOFC  anodes  is  also 
coupled  with  the  bulk  chemical  reactions  and  the  resulting  diffu¬ 
sion  and  reaction  time  scales  of  different  species  dictate  the  limiting 
behavior.  Accurate  modeling  of  this  phenomenon  is  imperative  for 
development  of  better  fuel  cell  designs. 

The  equations  governing  the  multi-component  mass  transfer 
include  Stefan-Maxwell  relations  for  calculation  of  specie  fluxes. 
These  relations  are  coupled,  nonlinear  system  of  partial  differen¬ 
tial  equations  and  are  cumbersome  to  solve,  especially  when  a 
large  number  of  species  are  involved  as  in  the  case  of  coal  syn¬ 
gas.  Alternatively,  simple  Fields  law  type  relations  with  an  effective 
multi-component  diffusivity,  which  can  be  directly  substituted  into 
the  specie  conservation  equation,  are  generally  used  to  model  the 
diffusion  fluxes  [1-4].  Flowever,  the  validity  of  such  simplified  mod¬ 
els  is  not  well  established.  Suwanwarangkul  et  al.  [5],  for  example, 
compared  these  two  approaches  for  a  simple  one-dimensional  case 
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without  bulk  chemical  reactions  and  recommended  the  use  of 
Stefan-Maxwell  relations  including  the  effect  of  Knudsen  diffusion 
for  multi-component  systems.  Flowever,  several  authors  [6-8]  pre¬ 
dicted  the  limiting  behavior  of  the  experimental  VI  curves  using 
Fields  law  type  relations  for  the  diffusion  fluxes.  In  such  studies, 
tortuosity  values  as  high  as  14.5  were  used  during  calibration  [6] 
whereas  some  studies  show  that  the  actual  tortuosities  for  porous 
electrodes  are  in  the  range  of  2-6  [9,10].  Thus,  it  is  necessary  to  look 
for  other  factors  that  may  affect  the  effective  multi-component  dif¬ 
fusion  coefficients  used  in  Fields  law  type  relations  rather  than  using 
unphysically  high  tortuosity  values. 

In  this  paper,  we  compare  these  two  approaches  by  simulat¬ 
ing  mass  transfer  inside  an  SOFC  anode  operating  on  coal  syngas 
consisting  of  CFI4,  H2,  CO,  FI2O,  C02  and  N2.  A  two-dimensional 
computer  code  was  developed  for  the  solution  of  models  based  on 
Stefan-Maxwell  and  Fields  law  flux  relations  to  predict  the  specie 
concentration  distribution  inside  the  porous  anode.  Methane  steam 
reforming  and  water  gas  shift  reactions  are  taken  into  account  in 
addition  to  the  electrochemical  reactions  at  the  anode/electrolyte 
interface  and  chemical  kinetics  modeling  was  employed  for  the  bulk 
reactions.  The  ratio  between  the  current  produced  by  CO  and  H2 
oxidation  is  assumed  to  be  constant.  Comparisons  show  that  at  low 
current  densities  the  models  agree  well  but  as  current  increases  the 
differences  increase.  The  error  in  the  Fields  law  model  was  positive 
for  the  reactants  and  negative  for  the  products.  Based  on  these  find¬ 
ings  an  empirical  correction  is  proposed  for  the  effective  diffusivity 
used  in  Fields  model.  Model  parameters  were  found  by  calibrating 
it  with  the  Stefan-Maxwell  model  using  the  two-dimensional  code. 
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Nomenclature 

c  model  parameter 

effective  diffusivity  (m2  s-1 ) 

Deff*  corrected  effective  diffusivity  (m2  s_1 ) 

Dj<  Knudsen  diffusivity  (m2  s_1 ) 

Dij  binary  diffusivity  (m2  s_1 ) 

£/<*  effective  Knudsen  diffusivity  (m2  s_1 ) 

Dt  effective  binary  diffusivity  (m2  s_1 ) 
i  total  current  density  (A  m- 2 ) 

iref  reference  current  density  (A  nrr  2 ) 

/<+  forward  reaction  rate  constant  for  methane  steam 
reforming  reaction  (molm-3  s-1  Pa-2) 
k~  backward  reaction  rate  constant  for  methane  steam 
reforming  reaction  (molm-3  s-1  Pa-4) 
kf,  kf  forward  and  backward  reaction  rate  constants  for 
water  gas  shift  reaction  (mol  m-3  s-1  Pa-2) 

Mi  molecular  weight  of  species  i  (g  mol-1 ) 

Mavg  average  molecular  weight  (g  mol-1 ) 
n  model  parameter 

N/  molar  flux  (mol  m-2  s-1 ) 

P  pressure  of  the  system  (Pa) 

R  universal  gas  constant  (J  mol-1  K-1 ) 

Ri  molar  production  rate  of  species  i  (mol  nrr3  s-1 ) 

(r)  mean  pore  radius  (pan) 

T  temperature  of  the  system  (K) 

j/j  mole  fractions 

s  porosity 

Vi  stoichiometric  coefficient 

r  tortuosity 


Proposed  correction  for  effective  diffusivity  was  then  applied  to  a 
three-dimensional,  in-house  SOFC  simulation  code  (DREAM-SOFC) 
which  uses  Fields  law  type  relations  for  diffusion  flux  calculations. 
As  an  application,  an  SOFC  button  cell  which  is  being  tested  at  West 
Virginia  University  was  simulated.  The  results  with  and  without  the 
proposed  correction  for  effective  diffusivity  were  compared. 


2.  Description  of  the  method 


Diffusion  transport  inside  the  porous  medium  has  two  com¬ 
ponents:  bulk  molecular  diffusion  where  molecular  to  molecular 
interaction  dominates  the  transport  and  Knudsen  diffusion  where 
molecular  to  solid  wall  interaction  dominates  the  transport.  Several 
different  approaches  have  been  developed  to  evaluate  multi- 
component  porous  media  transport  [3,11,12].  In  the  present  study, 
Stefan-Maxwell  model  (SMM)  and  Fields  model  (FM)  both  of  which 
include  the  effect  of  Knudsen  diffusion  in  addition  the  bulk  molec¬ 
ular  diffusion  were  used. 

The  flux  expression  in  Stefan-Maxwell  model  is  given  by  [13] 


Nj_  sr^yjNi  -yiNj 

[)k*  Dt 

j±i 


(i) 


where  N,-  and  y2-  are  the  molar  fluxes  and  mole  fractions  of  species, 
Dt  and  are  effective  binary  and  Knudsen  diffusivities.  The  first 
term  on  the  left  hand  side  accounts  for  the  Knudsen  diffusion  and 
the  second  one  signifies  the  multi-component  bulk  molecular  dif¬ 
fusion. 


The  effective  binary  diffusivity  is  usually  defined  as: 


Dt  =  -Du 


(2) 


where  £  and  r  are  the  bulk  porosity  and  tortuosity,  respectively.  Dy  is 
the  binary  diffusivity  and  can  be  calculated  using  Chapman-Enskog 
relations  [13]: 

,1/2 


Dy  =  0.001858 


[T3(M,-  +  Mj)/MjMj] 
PcrfjQD 


(3) 


where  T  and  P  are  the  temperature  and  the  pressure  of  the  system, 
respectively.  M,-  represents  the  molecular  weight  of  species  i.  cry  is 
given  by 


Gi+Gj 


(4) 


a/  is  the  collision  diameter  of  species  i.  The  dimensionless  param 
eter  £2D  is  the  collision  integral  calculated  as: 

CFG 


o  _  A 

D  ~  T%  +  exp (DTW) 


exp(FTN)  exp  (HTn) 


(5) 


where  the  coefficients  A,  B,  C,  D,  E,  F,  C,  H  are  1.06036,  0.1561,  0.193, 
0.47635,  1.03587,  1.52996,  1.76474,  3.89411,  respectively  [14].  The 
dimensionless  temperature  TN  is  defined  as: 


blJ 

where  k  is  the  Boltzmann  constant.  The  term  £y  is  calculated  as: 


(6) 


i  =  ( SiSj ) 


1/2 


(7) 


where  £z  is  the  characteristic  Lennard-Jones  energy  of  species  i. 
Likewise  Knudsen  diffusivity  is  defined  as: 


(8) 


where  (r)  is  the  mean  pore  radius,  R  is  the  universal  gas  constant 
and  Mi  is  the  molecular  weight  of  the  species  i. 

Fields  model  (FM),  simplest  form  used  to  describe  the  transport 
of  species  in  the  gas  phase,  is  defined  as: 

N,  =  -Df  F-Vy,  (9) 

where  D*f  is  the  effective  diffusivity.  Simulations  with  Fields  model 
were  conducted  by  using  the  effective  diffusivity  proposed  by  Yak- 
abe  et  al.  [3]: 


jjeff  _  £  1  Wj^yj  J_ 

1  -  r  [  Di  m  +  D* 

ofj?m  is  a  dimensionless  parameter  given  by 


(10) 


—  1 


Mi 

Mavg 


1/2 


(11) 


where  Mavg  is  the  average  molecular  weight  of  the  mixture  and  Mi  is 
the  molecular  weight  of  the  species  i.  Multi-component  molecular 
diffusivity  of  species  i  in  the  mixture  Dim  is  given  by 

DUm  =  (12) 

*  i  Dy 

The  calculated  molar  fluxes  either  from  SMM  (Eq.  (1))  or  from  FM 
(Eq.  (9))  are  then  used  to  determine  the  rate  of  change  of  specie 
concentrations  through  the  following  conservation  equation: 


ptyt 

RT  dt 


=  -VNj  +  Rj 


(13) 
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x  =  0  x  =  0.542  mm 

Fig.  1.  Schematic  representation  of  components  of  a  solid  oxide  fuel  cell. 


where  Rf  is  the  molar  rate  of  production  of  species  i  due  to  bulk 
chemical  reactions. 

In  the  present  study,  the  principle  bulk  chemical  reactions  of 
concern  are  methane  steam  reforming  and  water  gas  shift  reactions, 
described  respectively  as: 

k+ 

CH4  +  H20^3H2  +  CO  (14) 

kr 

C0  +  H20^C02  +  H2  (15) 

K 

where  kf  and  kf  are  the  forward  and  backward  reaction  rate  con¬ 
stants  for  reforming  reaction  and  /<+  and  kf  are  the  reaction  rate 
constants  for  shift  reaction.  Experimentally  determined  values  for 
these  rate  constants  and  their  Arrhenius  form  can  be  found  in 
[15,16].  The  net  reaction  rates  are  given  for  reforming  and  shift 
reactions,  respectively  as: 

Rr  =  /<^Pch4Ph2o  -  /<r“Pco(PH2)3  (16) 

Rs  =  kfpcoPu2o  ~  KPco2Ph2  (17) 

The  specified  boundary  conditions  are  the  specie  concentra¬ 
tion  at  the  anode  fuel  channel  interface  and  molar  fluxes  at 
anode-electrolyte  interface  proportional  to  the  prescribed  cell  cur¬ 
rent.  It  was  assumed  that  there  is  zero  flux  for  other  boundaries. 
The  ratio  between  the  current  produced  by  CO  and  H2  oxidation  is 
assumed  to  be  constant  and  taken  as  4  for  the  current  study  [17]. 

3.  Results  and  discussions 

As  explained  in  the  previous  section  SMM  and  FM  are  two 
of  the  models  which  are  used  to  calculate  the  fluxes  of  each 
species  in  a  multi-component  system.  In  this  paper,  these  two 
approaches  were  compared  by  simulating  mass  transfer  inside  an 
SOFC  anode  operating  on  coal  syngas  consisting  of  CH4,  H2,  CO, 
H20,  C02  and  N2  (see  Table  3).  A  schematic  representation  of  the 
SOFC  anode  which  is  used  for  simulation  is  shown  in  Fig.  1.  A 
two-dimensional  computer  code  was  developed  for  the  solution 
of  models  based  on  Stefan-Maxwell  and  Fields  law  flux  relations 
to  predict  the  specie  concentration  distribution  inside  the  porous 
anode. 

The  predictive  ability  of  the  developed  code  is  verified  by  com¬ 
paring  the  results  of  molar  formation  rates  to  those  of  [  1 5  ]  that  also 
used  SMM  to  predict  the  species  molar  fluxes.  Fig.  2  shows  a  good 
agreement  between  two  independent  calculations  which  indicates 
the  model  was  successfully  implemented. 

In  what  follows,  the  2D  code  validated  above  is  used  to  compare 
the  results  obtained  using  Fields  model  for  diffusion  flux  with  those 
obtained  using  Stefan-Maxwell  model.  Fig.  3  shows  the  maximum 


Electrode  Position  (m) 

Fig.  2.  Comparison  between  the  molar  rate  of  formation  results  of  the  present  study 
and  those  of  [15]. 


percentage  error  (deviation  from  SMM  predictions  of  concentra¬ 
tion)  of  Fields  model  versus  current  density. 

As  can  be  seen  from  the  figure,  at  low  current  densities  the  mod¬ 
els  agree  well  but  as  current  increases  the  differences  increase  to 
up  to  20%  at  a  current  density  of  4  A  cm-2.  The  error  in  the  Fields 
law  model  is  positive  for  the  fuel  species  (reactants)  and  negative 
for  the  products.  Based  on  these  findings  an  empirical  correction  is 


Fig.  3.  Maximum  percentage  error  of  Fields  model  versus  current  density  for  (a) 
reactants,  (b)  products. 
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Table  1  Table  2 

Coal  syngas  composition.  Anode  properties. 


Species 

Mole  fractions 

Value 

CO 

0.291 

Thickness  (mm) 

0.5 

h2 

0.285 

Width  (mm) 

0.542 

co2 

0.118 

Tortuosity 

3.6 

h2o 

0.276 

Porosity 

0.56 

ch4 

0.021 

Mean  pore  diameter,  (r)  (p,m) 

1.07 

n2 

0.009 

Operating  temperature  (K) 

1123 

proposed  for  the  effective  diffusivity  used  in  Fields  law  model  as: 

Dff*  =  Df  ^1  +sign(uI)c^-^  (18) 

where  v*  is  the  stoichiometric  coefficient,  i  is  the  total  current  den¬ 
sity  and  iref  is  the  reference  current  density  taken  as  1  A  m-2.  Model 
parameters  c  and  n  are  found  as  4.88  x  10-4  and  2,  respectively, 
by  calibrating  it  with  the  SMM  using  the  two-dimensional  code 
presented  here. 

The  corrected  diffusivities  were  used  to  predict  the  species 
concentration  distribution  inside  an  SOFC  anode  running  on 
coal  syngas.  Anode  was  treated  as  a  two-dimensional  medium. 
Tables  1  and  2  show  the  syngas  composition  and  the  anode  elec¬ 
trode  properties  used  for  simulations. 

Mole  fraction  distributions  inside  the  porous  anode  directly 
below  the  fuel  channel  is  given  in  Figs.  4  and  5  for  operating  current 
densities  of  1 A  cm-2  and  4  A  cm-2 ,  respectively.  As  can  be  seen  from 
the  figures  the  corrected  diffusivity  coefficient  gave  better  agree¬ 
ment  with  the  Stefan-Maxwell  model  and  even  at  higher  current 


Table  3 

Normalized  diffusion  time  scales. 


densities  the  error  is  less  than  5%  (see  Fig.  3).  Since  this  correction  is 
based  on  numerical  experiments  it  is  expected  to  have  some  errors 
and  5%  error  is  acceptable  for  such  a  heuristic  model. 

The  estimated  reaction  and  diffusion  time  scales  for  the  con¬ 
ditions  used  in  the  above  simulation  normalized  with  the  flow 
through  time  (r Conv  =  Llu  =  0.267  s  with  L  being  fuel  channel  length 
and  u  being  average  velocity  inside  the  channel)  of  the  fuel  in  the 
gas  channel  are  given  in  Tables  3  and  4  respectively.  It  can  be  seen 
from  Table  3  that  there  is  ample  time  for  the  component  species 
to  diffuse  through  the  anode  as  they  flow  through  the  gas  channel. 


Fig.  4.  Concentration  profiles  along  the  anode  thickness  predicted  by  FM  and  SMM  for  (a)  CO,  (b)  H2,  (c)  C02  and  (d)  H20  at  1 A  cm-2. 
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Fig.  5.  Concentration  profiles  along  the  anode  thickness  predicted  by  FM  and  SMM  for  (a)  CO,  (b)  H2,  (c)  C02  and  (d)  H20  at  4  A  cm-2. 


Normalized  reaction  time  scales  (Table  4)  indicate  that,  with  the 
given  inlet  concentrations,  steam  methane  reformation  mostly  pro¬ 
ceeds  in  the  forward  direction  whereas  as  water  gas  shift  reactions 
may  proceed  in  both  directions.  Moreover  diffusion  time  scales  are 
much  smaller  compared  to  reaction  time  scales  indicating  that  this 
electrode  is  not  limited  by  diffusion. 

4.  Application:  button  cell 

In  order  to  independently  verify  the  proposed  corrections  to 
the  effective  diffusivities  used  in  Fields  law  type  models  applied 
to  the  multi-component  mass  transfer  problems  in  fuel  cells;  this 
approach  is  employed  in  a  three-dimensional  simulation  of  a  real 
button  cell.  The  model  used  for  these  simulations  is  described  else¬ 
where  [18,19]  in  the  literature  and  will  not  be  covered  here.  The 
salient  features  of  the  model  are  that  it  allows  for  several  species 
and  bulk  reactions  on  the  anode  side  and  it  handles  the  simulta¬ 
neous  electrochemical  oxidation  on  CO  and  H2  without  assuming 
a  constant  ratio  between  the  current  produced  by  each;  rather  this 
ratio  is  an  outcome  of  the  model  (see  Appendix).  Also,  the  diffu¬ 
sion  fluxes  are  calculated  using  Fields  law  type  relations  and  the 
effective  diffusivities  are  calculated  using  the  model  in  Eq.  (10).  The 


Table  4 

Normalized  reaction  time  scales. 


Tchemltconv 

Forward 

Backward 

Steam  methane  reforming 

0.169621 

1570.154 

Water  gas  shift 

0.069015 

0.048962 

button  cell  simulated  here  is  modeled  after  the  experiments  being 
performed  at  WVU  [20].  The  active  area  of  the  button  cell  is  approx¬ 
imately  2  cm2.  The  anode,  electrolyte  and  cathode  thicknesses  are 
780  |jim,  20  |jim  and  200  p,m  respectively.  The  fuel  stream  fed  to  the 
anode  has  the  composition  shown  in  Table  5.  The  cathode  side  is 
supplied  with  air  (21%  02  and  79%  N2).  Both  air  and  fuel  streams 
are  supplied  at  a  temperature  of  800  °C.  The  flow  field  around  the 
button  cell  is  not  included  in  the  computation  domain.  However, 
its  effect  on  the  cell  is  taken  into  account  using  appropriate  bound¬ 
ary  conditions  and  transport  properties.  Given  the  low  utilization 
of  fuel  and  air  in  the  button  cell  experiments,  it  was  found  that 
the  current  density  distribution  inside  the  cell  was  uniform  at  all 
currents. 

Fig.  6  shows  the  experimentally  measured  and  numerically  com¬ 
puted  VI  curves  for  the  button  cell.  Computations  are  performed 
both  with  and  without  the  proposed  correction  to  the  effective 
multi-component  diffusion  coefficients.  It  can  be  seen  from  Fig.  6 
that  the  agreement  between  the  experiments  and  numerical  pre¬ 
dictions  is  good  for  low  cell  currents  but  at  high  currents  the 
computations  without  correction  to  the  diffusion  coefficients  fail 
to  predict  the  onset  of  diffusion  limitation  of  the  cell  performance. 


Table  5 

Fuel  composition  supplied  to  the  anode. 


Species 

Mole  fractions 

h2 

0.303 

CO 

0.229 

co2 

0.209 

H20 

0.259 
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Fig.  6.  Experimental  and  numerically  simulated  VI  curves  for  the  button  cell. 


However  the  numerical  results  with  corrected  diffusion  coefficients 
accurately  predict  the  limiting  current  though  the  agreement  in 
terms  of  voltage  values  is  poor  near  the  limiting  current  region  of 
the  VI  curve.  It  should  be  noted  that  a  reasonable  value  of  2.5  [9,10] 
is  used  for  tortuosity  in  these  simulations.  Though  some  researches 
[21,22]  use  higher  values  for  tortuosity  to  calibrate  diffusion  limi¬ 
tation,  they  are  not  believed  to  be  physically  realistic  [9,23].  For  a 
thorough  analysis,  the  simulations  were  repeated  with  a  higher  tor¬ 
tuosity  of  10  with  and  without  correction  and  the  resulting  VI  curves 
are  also  plotted  in  Fig.  6.  It  can  be  seen  that,  although  the  voltages 
for  r  =  10  are  lower  than  those  for  r  =  2.5,  the  curve  does  not  show 
the  typical  drop  in  voltage  at  high  currents  when  no  correction  was 
applied.  When  the  correction  was  applied  with  r  =  10,  the  cell  lim¬ 
iting  current  was  found  to  be  much  lower  than  the  experimentally 
observed  value.  It  may  thus  be  concluded  that  the  inability  of  most 
numerical  models  to  accurately  predict  the  diffusion  limitation  is 
due  to  the  inaccuracy  resulting  from  Fields  law  type  diffusion  model 
rather  than  the  uncertainty  in  the  tortuosity  values.  The  limiting 
behavior  is  known  to  be  due  to  the  concentration  of  the  reactants 
becoming  very  small  near  the  active  electrolyte/electrode  interface. 

In  order  to  see  how  the  concentrations  of  the  H2  and  CO  near  the 
active  interface  are  varying  with  the  cell  current,  their  mole  frac¬ 
tions  are  plotted  against  the  average  cell  current  density  in  Fig.  7. 
Here  it  is  seen  that  there  is  a  huge  difference  in  the  prediction 
of  interface  concentrations  obtained  from  simulations  with  and 
without  the  proposed  correction,  especially  at  high  currents.  The 
interface  concentrations  in  the  case  of  simulations  without  correc¬ 
tion  hardly  change  with  cell  current  when  compared  to  those  with 
the  correction.  This  is  due  to  the  reduction  in  diffusion  coefficient 
as  a  function  of  current  through  the  proposed  correction.  This  could 
have  significant  effect  on  the  concentration  distributions  and  also 
chemical  kinetics  inside  the  anode.  For  example,  the  profiles  of  H2 
and  CO  mole  fractions  along  the  normalized  thickness  of  the  anode 
are  shown  in  Figs.  8  and  9  respectively  for  the  case  of  average  cell 
current  density  1 A  cm-2.  Also  the  local  production  rate  of  H2  (due 
to  bulk  reactions  inside  the  anode)  along  the  normalized  thickness 
of  anode  for  the  same  case  is  shown  in  Fig.  10.  It  can  be  seen  that  all 
these  profiles  are  significantly  different  for  cases  with  and  without 
the  correction.  More  interestingly,  the  production  of  H2  is  nega¬ 
tive  for  the  case  without  correction,  indicating  net  consumption  of 
hydrogen  in  the  water  gas  shift  reaction  whereas  it  is  positive  for 
the  case  with  correction.  These  differences  are  less  prominent  for 
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Fig.  7.  Concentrations  of  H2  and  CO  at  the  active  anode/electrolyte  interface  as  a 
function  of  cell  current  density. 

cell  currents  lower  than  the  case  shown  here  and  more  prominent 
for  higher  cell  current.  These  differences  are  to  be  expected  since 
the  concentrations  and  reaction  rates  are  interdependent  on  each 
other  and  are  sensitive  to  the  diffusion  coefficient  which  is  widely 
different  between  these  two  cases,  especially  at  high  currents. 

Fig.  11  shows  the  ratio  of  the  current  produced  by  H2  to  the  cur¬ 
rent  produced  by  CO  oxidation  as  a  function  of  cell  current  density. 
As  was  mentioned  before  this  ratio  is  not  assumed  to  be  constant  as 
done  by  some  authors  [17,24,25]  but  was  determined  as  part  of  the 
solution.  It  can  be  seen  from  Fig.  11  that  this  ratio  is  in  fact  a  func¬ 
tion  of  the  cell  current  density  and  as  with  the  other  variables,  its 
value  is  affected  by  the  proposed  correction.  With  corrected  diffu¬ 
sion,  the  proportion  of  current  from  CO  oxidation  seems  to  decrease 
with  the  cell  current  which  is  a  result  of  nearly  zero  concentration 
of  CO  near  the  active  interface  at  these  currents  as  can  be  seen 
from  Fig.  7.  It  should  also  be  noted  that  in  the  simple  mass  transfer 


Fig.  8.  Concentrations  of  H2  along  the  normalized  thickness  of  the  anode  for  cell 
current  density  of  1 A  cm-2. 
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Fig.  9.  Concentrations  of  CO  along  the  normalized  thickness  of  the  anode  for  cell 
current  density  of  1 A  cm-2. 

calculations  used  to  formulate  the  proposed  correction,  a  constant 
ratio  of  4  is  assumed  between  H2  and  CO  currents  whereas  this  ratio 
seems  to  be  getting  as  high  as  7  when  this  correction  is  applied  in  a 
real  cell.  It  remains  to  be  tested  if  the  correction  is  still  valid  under 
such  conditions  where  this  ratio  is  a  function  of  cell  current  as  well. 
This  could  be  a  possible  reason  for  the  disagreement  between  the 
voltage  predictions  using  the  corrected  diffusion  coefficient  and  the 
experiments  shown  in  Fig.  6. 

Finally,  to  investigate  the  effect  of  different  fuel  compositions 
on  the  limiting  current,  button  cell  simulations  are  repeated  using 
r  =  2.5  with  corrected  effected  diffusivities  and  the  results  are 
shown  in  Fig.  12.  The  model  results  show  expected  dependence  of 
VI  curves  on  the  fuel  concentrations  with  voltages  and  limiting  cur¬ 
rents  being  higher  when  the  concentration  of  fuel  species  (CO  and 
H2)  is  higher.  Thus,  it  is  shown  that  the  proposed  correction  does 
not  alter  the  response  of  the  model  to  changes  in  fuel  composition. 


Fig.  10.  Local  production  rate  of  H2  along  the  thickness  of  the  anode  for  cell  current 
density  of  1 A  cm-2. 


Fig.  11.  Ratio  between  the  current  produced  from  H2  oxidation  and  that  from  CO 
oxidation  as  a  function  of  cell  current. 


5.  Conclusions 

Stefan-Maxwell  model  (SMM)  and  simple  Fields  model  (FM) 
type  of  relations  both  including  Knudsen  diffusion  for  the  cal¬ 
culation  of  species  mole  fractions  in  a  multi-component  systems 
were  compared  and  it  was  found  that  at  low  current  densities  the 
models  agree  well  but  as  current  increases  the  differences  also 
increase.  Based  on  the  findings  an  empirical  correction  is  proposed 
for  the  effective  diffusivity  used  in  Fields  model.  The  corrected  diffu- 
sivity  coefficient  gave  better  agreement  with  the  Stefan-Maxwell 
model  and  even  at  higher  current  densities  the  error  is  less  than 
5%.  Proposed  correction  for  effective  diffusivity  was  applied  to 
a  three-dimensional,  in-house  SOFC  simulation  code  which  uses 
Fields  model  type  relations  for  diffusion  flux  calculations.  As  an 
application,  an  SOFC  button  cell  which  is  being  tested  at  West  Vir¬ 
ginia  University  was  simulated.  Numerical  simulations  with  the 
corrected  diffusivity  accurately  predicted  the  limiting  current  even 
with  a  reasonable  tortuosity  of  2.5.  However,  the  predicted  voltages 
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were  slightly  different  from  the  experimental  values  near  the  lim¬ 
iting  current  region.  Also  the  distributions  of  concentrations  and 
reaction  rates  are  widely  different  for  the  predictions  with  and 
without  the  proposed  correction,  especially  at  high  current  density. 
The  ratio  between  the  currents  produced  from  H2  and  CO  oxida¬ 
tion  is  a  function  of  cell  current  density  and  is  also  affected  by  the 
correction. 
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Appendix  A 


Electrochemistry  model 

The  half  cell  reactions  considered  in  the  anode  are 


H2  +  CT^  H20  +  2e- 

(Al) 

CO  +  0=^  C02  +  2e- 

and  in  the  cathode  the  half  cell  reaction  is 

(A2) 

02  +4e_  2CT 

Two  overall  reactions  can  be  established  as: 

(A3) 

H2  + 1  /202  H20 

(A4) 

CO  +  1  /202  C02 

(A5) 

The  Nernst  potential  for  reactions  ( A4)  and  ( A5 )  is  given  respectively 

by 


£Hp  =  £„  -  ^  ln 
H?  h2  2  F 


q  RT  . 
Em  —  Eq0  -  —  In 


2  F 


Ph2o 

Ph2PS25 

Pco2 

PcoPf°25 


(A6) 


(A7) 


Following  the  idea  proposed  by  Nishino  et  al.  [26]  the  potential 
across  the  electrolyte  is  given  by 


Eh2  ~  hn2  ~  ho2  =  Ec  -  Ea 
Em  ~  dco  ~  ho2  =  Ec  -  Ea 


(A8) 

(A9) 


The  calculation  of  the  overpotentials  (z7h2,  hco  and  r]o2)  is  per¬ 
formed  by  using  the  Butler-Volmer  equation: 


h  =  io,s 


exp 


fansF 


%)  -  exp  ( 


(1  -a)nsF 


V  RT  'ISJ  6XP  V  RT  >h 


(A10) 


where  the  subscript  s  represents  the  species  H2,  CO,  and  02. 


itot  =  i  H2+ico  (All) 

Eqs.  (A8)-(A11 )  form  a  set  of  equations  to  solve  for  hydrogen  and 
carbon  monoxide  activation  overpotentials  and  current  densities. 
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